Impact ionization fronts in Si diodes: Numerical evidence of superfast propagation 

due to nonlocalized preionization 
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We present numerical evidence of a novel propagation mode for superfast impact ionization fronts 
in high-voltage Si p+-n-n+ structures. In nonlinear dynamics terms, this mode corresponds to a 
pulled front propagating into an unstable state in the regime of nonlocalized initial conditions. Before 
the front starts to travel, field-ehanced emission of electrons from deep-level impurities preionizes 
initially depleted n base creating spatially nonuniform free carriers profile. Impact ionization takes 
place in the whole high-field region. We find two ionizing fronts that propagate in opposite directions 
with velocities up to 10 times higher than the saturated drift velocity. 
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Superfast impact ionization fronts travel with veloc- 
ity vf higher than the saturated drift velocity Vg. Ex- 
citation of such fronts is the fastest non-optical method 
to modulate conductivity of a high- voltage semiconduc- 
tor structurei^"— This method has important pulse-power 
applicationi^""— Recently, we proposed a novel mode of su- 
perfast front propagation in semiconductor structures: a 
pulled front propagating into unstable state in the regime 
of nonlocalized initial conditions^ This mode is expected 
to appear in p'^-n-n'^ structures with a relatively low n 
base doping level. The propagation mechanism is qualita- 
tively different from that for the well-known TRAPATT- 
like propagation mode^iSr— widely considered as a most 
feasible mechanism of the superfast switching of high- 
voltage pulse-power devicesj^ For the impact ionization 
front the term "nonlocalized initial conditions"'iS. means 
a small free carriers concentration that decreases in the 
direction of front propagation. It has been suggested^ 
that such preionization of the initially depleted n base 
may be created by field-enhanced electron emission from 
deep-level centers.— Numerical simulations presented in 
the Letter confirms this suggestion and provide the first 
evidence of superfast front propagation due to nonlocal- 
ized preionization in realistic Si p^-n-ri^ structures used 
in pulse-power applications. 

Process-induced (PI) deep-level centers in Si are 
double-charged donors with ionization energies 0.28 eV 
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FIG. 1: Sketch of the reversely biased p 
erated in the external circuit. 
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of as recombination centers,— PI centers turned out to 
be deep electron traps. ^ As they do not influence the 
life-time of nonequilibrium carriers (the most carefully 
controlled parameter of the commercial material), their 
presence in high-purity Si is "hidden" .— In high voltage 
structures used in power applictations PI centers appear 
in concentration iVpi = 10"... lO^^ cm-^ (Ref.[ll). Field- 
enhanced ionization of PI centers is a potential mecha- 
nism of determenistic low-jitter triggering of the ionizing 
front in high-voltage structures.^ In the simplest case of 
low temperatures (T < 200 K) the emission rate depends 
on the electrical field F as(Ref. 
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where Fq = Ay/2mE^/3qh, Eq ^ 0.28 eV is the binding 
energy (U level), Eb is Bohr energy in semiconductor, 
m is the effective electron mass and q is the elementary 
charge. Field-enhanced ionization of PI centers has a 
characteristic threshold F^^^ w 3 • 10^ V/cm (Ref'?) that 
exceeds the threshold of band-to-band impact ionization 
F™P « 2 • 10^ V/cm. It is important that F^^^ > F™p 
because electrical field in the n base should be increased 
above J^™^ before impact ionization starts,^ 

Superfast triggering occurs when a sharp voltage ramp 
is applied to reversely biased Si p^-n-n"*" structure con- 
nected to the voltage source V{t) in series with load re- 
sistance R (Fig. 1). The applied voltage is modelled as 
V{t) = Vq + A ■ t, where the inital voltage Vq is chosen 
so that the n base is fully depleted at t = 0. We assume 
that p'^-n and n-n"*" junctions are sharp and restict mod- 
elling to the n base. To describe the dynamics of elec- 
tron and hole concentrations n{z, t), p{z, t) and electrical 
field F{z,t), we use the standard one-dimensional drift- 
diffusion model together with the Poisson equation and 
the Kirchoff equation for the external circuit. The gen- 
eration term includes the cut-off that eliminates unphys- 
ical multiplication at low concentrations n,p < ricut^ 
We refer to Ref. [l^ for the details of the model and of 
the numerical method used. The generation term addi- 




FIG. 2: Voltage U(t) during the switching process for differ- 
ent voltage ramps A = 0.5, 1, 5 kV/ns (curves 1,2 and 3, 
respectively) . 



tionally incorporates field-enhanced emission of free elec- 
trons from PI centers with a rate given by Eq. ([!]). We 
use the follovifing set of structure and circuit parameters: 
the n base length W = 100 /im, dopant concentration 
iVd = lO^^cm"^, A'pi = lO^^cm"'^, cross-section area 
S = 0.02 cm~^, load resistance R = 50 Ohm, initial bias 
Vq = 1 kV. The respective stationary breakdown voltage 
is about 1.2 kV. We choose ricut = 10^ cm~^, as in Refs. 
[llIlO. 

In Fig. 2 we show the voltage U (t) over the diode dur- 
ing switching process. At the first stage U {t) nearly fol- 
lows the applied voltage V{t). Then, U{t) sharply drops 
and the current flowing through the diode increases. The 
triggering time is close to 100 ps and decreases with in- 
crease of the voltage ramp A. This is 10 times faster than 
the drift time W/vg « 1 ns. Such superfast switching has 
been observed in the whole range of actual voltage ramps 
0.5. ..10 kV/ns. Although the transient U{t) looks similar 
to the earlier numerical results for TRAPATT-like^ii^-iS 
or tunneling-assisted^Sii^ impact ionization fronts, the in- 
ner dynimics turns out to be qualitatively different. 

In Fig. 3 we show the inner dynamics for A = 1 kV/ns. 
Due to the higher electrical field, the emission of electrons 
from PI centers is most efficient near p'^-n junction. Free 
electrons drift to the left and multiplicate by the band-to- 
band impact ionization. Due to this drift, the maximum 
of the concentration profile is shifted from the p~^-n junc- 
tion into the n base [Fig. 3(a)]. Screening begins when 
concentrations n and p overcome N^. Eventually, the 
avalance multiplication creates the initial plasma layer 
that fully screens the applied electrical field [Fig. 3(a), 
curve 3]. Afterwards, the plasma layer expands due to 
propagation of two ionizing fronts travelling in opposite 
directions [Fig. 3(b)]. 

The left and the right (negative and positive) fronts 
travel with velocities vj « 10 Vs and v'^ « 3Vs, respec- 
tively. Both fronts propagate into the areas where the 
electrical field i^m is nearly constant in space and varies 
in the interval from 3 • 10^ to 4 • 10^ V/cm, nearly 2 times 
larger than i^™'^- The small amount of free cariers is 
present everywhere (Fig. 4). Hence avalance multiplica- 
tion also goes on all over the high-field region. These con- 





FIG. 3: The spatial profiles of the electrical field F{z, t) 
and electron and hole concentrations n{z,t) and p(z,t) in 
the n-base at different times: (a) increase of the free car- 
rier concentration due to the field-enhanced ionization of 
deep-level impurities and subsequent avalanche multiplica- 
tion at times t = 2.285, 2.295, 2.31 ns (curves 1,2, and 3, re- 
spectively); (b) propagation of impact ionization fronts at 
t = 2.315, 2.325, 2.335, 2.345, 2.355 ns (curves 1,2,3,4, and 5, 
repectively). Numerical parameters as in Fig. 2 for A = 
1 kV/ns. 




FIG. 4: The spatial profile of electron concentrations n{z, t) 
in the propagating front at times t — 2.315, 2.325, 2.335 ns 
(curves 1,2, and 3, respectively,) in logarithmic scale. Numer- 
ical parameters as in Fig. 3. 



ditions are favorable for quasiuniform breakdown. The 
reason why the front-like solution nevertheless occurs is 
the nonuniform profile of preionization: the concentra- 
tion of free carriers in the high-field region decreases in 
the direction of front propagation. This decay happens 
to be nearly exponential n,p ~ exp(±Az) with character- 
istic exponents A~ « 4-10'^ cm~^ and A"*" w 1.3-10^ cm~^ 
for negative and positive fronts, respectively (Fig. 4). 
In nonlinear dynamics, such fronts are known as pulled 
fronts propagating into unstable state in the regime of 
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nonlocalized initial conditions'^ 

Let us compare the numerical results with analytical 
predictions for front velocity and plasma concentration 
made in Ref. for the simplified model where electrons 
and holes have the same impact ionization coefficients 
a{F) = an{F) = ap{F) and Nd = 0. For sufficiently 
small exponent A, the front velocity is given by 

Vf/vs w 2a,„/A, (2) 

where am = C({F^) is the impact ionization coefficient 
in the high field region the front propagates to. In Si 
in strong electrical fields impact ionization by holes is 
much weaker than by electrons. Therefore 2ani shall be 
replaced by a„i = an(-Fin) in the nominator of Eq. ([2|). 
The numerical value is am ~ 3.2 • 10'' cm~' for Fm = 
3.5 • 10^ V/cm. Eq. © is apphcable if A <C A*, where A^ 
corresponds to so called marginally stable frontH and is 
explicitly determined by am, diffusion coefficient D and 
Vs (see Ref. d). Simple calculation shows that A* w 
2 • 10^ cm~^. Hence, the condition A* ^ A* is met for 
the case under consideration, and we obtain analytical 
estimates vj w am/A~ 8vs and v'j' sa am/A"*" sa 2.5 
for the negative and positive fronts shown in Figs. 3,4, 
respectively, in good agreement with numerical results. 

For the symmetrical model, the plasma concentra- 
tion a = n + p behind the front is given by tTpi = 

(^seQaQFi^l q) " " exp(— l/a;)da;, where the ionization 
coefficient is approximeted as a{F) — aoexp(—Fo/F) 
(Ref. 8). Removing again the factor of 2 and taking the 
numerical values ao — 7.4 • 10^ cm, Fo = 1.1 • 10^ V/cm 



for electrons in Si, we get (Jpi ~ 10'^ cm~^ for the actual 
electrical field Fm = 3.5 • 10^ V/cm, in good agreement 
with numerical results. 

The velocity of the convential TRAPATT-like front is 
determined by the size of impact ionization zonciiSr— In 
contrast, the velocity of the pulled front in the regime of 
nonlocalized initial conditions is determined by the slope 
of the preionization profile, while the size of impact ion- 
ization zone coincides with the total size of the high-field 
region. Consequently, the pulled fronts shall be expected 
in structures with moderate doping Nd where the condi- 
tion F > is met all over the high-field region. How- 
ever, a certain level of the n base doping Nd is crucial 
because it leads to the slope of the electrical field in the 
depleted n base. Due to this slope the field-enhanced ion- 
ization of PI centers is spatially nonuniform and creates 
decaying profile of preionization. In particluar, in p-i- 
n structures nonlocalized initial conditons for the pulled 
front propagation cannot be created by deep-level im- 
purities. For large Nd the convential TRAPATT-like 
modoiSi"— occurs, but the cross-over to the pulled mode 
is possible at the last stage of front propagation when the 
condition F > -FJ.™^ is met everywhere. 
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